The structure of the strange sea of the proton 
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In this work we study the strange sea of the proton using a version of the Meson Cloud Model 
containing both, efective and perturbative degrees of freedom. We construct the s and s parton 
distributions functions at the initial energy scale, Qq, where QCD evolution starts. The initial s and 
s pdfs depend on a number of parameters which we fix by comparison to parameterizations of the 
strange sea of the nucleon obtained in a recent global fit to experimental data, allowing for a s — s 
asymmetry. We show that the model describes well the strange sea of the proton and argue that 
it can be a phenomenologically motivated alternative to the usual input parameterizations used in 
fits to experimental DIS data. 
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I. INTRODUCTION 

The composite nature of hadrons in term of quark and 
gluon degrees of freedom is today firmly stablished. It is 
also firmly believed that the internal dynamics of hadrons 
is determined by the strong interactions between quarks 
and gluons, as governed by Quantum Chromodynamics 
(QCD). However, a detailed theoretical description of 
the hadron structure is still missing because QCD can 
only be solved in the perturbative regime, correspond- 
ing to the short distance domain probed in hard colli- 
sions, whereas the long distance part of the interaction 
requires a non-perturbative treatment usualy supplied by 
effective models, latice simulations, etc. It is worth not- 
ing that is just the long distance part of the interaction 
which is responsible for the hadron as a bound state of 
quarks and gluons. Indeed, hadrons are made up of a 
fixed number of valence quarks plus a varying number 
of sea quark-antiquark pairs and gluons, being these last 
the "glue" which keeps valence quarks together forming 
a hadron. As a matter of fact, the sea q — q pairs and 
gluons which bring together the valence quarks into a 
bound state are the so-called intrinsic sea [T| , which has 
to be distinguished from the extrinsic sea generated by 
QCD evolution and consequently, dependent on the en- 
ergy scale Q^. 

From an experimental point of view, the struc- 
ture of nucleons is by far the best known hadron 
structure, existing today a variety of parton distribu- 
tion function (pdf) parameterizations extracted from 
data. Although the analysis of the existing data has 
confirmed the impresive success of the Quark Par- 
ton Model and the Dokshitzer-Grivov-Lipatov-Altarelli- 
Parisi (DGLAP) evolution equations [2] when describing 
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the nucleon structure, in the Bjorken regime, in terms of 
pdfs given as a function of x, the fraction of the nucleon 
momentum carried by partons, and , the energy scale 
which in Deep Inelastic Scattering (DIS) experiments is 
identified with the 4-momentum transfer between the lep- 
ton and the nucleon, several questions remain still un- 
solved. Among them, what is the functional form - and 
why - of the valence and sea quark and gluon pdfs at the 
energy scale Qq where QCD evolution starts; what is the 
value of the initial Qp energy scale; what is the structure 
of the so-called intrinsic sea oi q — q pairs and gluons; etc. 
Of course, finding answers to these questions requires to 
deal with the long distance - confining - realm of QCD, 
thus, as long as there is no solution of QCD in the low 
regime, we have to rely on effective models. 

Furthermore, the structure of the nucleon's intrinsic 
sea of quarks and gluons reflects the dynamics of non- 
perturbative QCD. In fact, the form of the initial pdfs 
at the Qq initial scale, as well as the relative abundance 
of the different quark flavors are the footprint of sub- 
tle non-perturbative QCD dynamical effects. Notice, for 
instance, that the experimentally observed d/u asym- 
metry ^ and the violation of the Gottfried sum rule 
(GSR) 4] cannot be described in terms only of perturba- 
tive QCD and gluon splitting [S], as the ratio d/u would 
became equal to one, due to the equal probability of 
gluon splitting into dd or uu pairs. The distribution of 
strange and anti-strange quarks in the nucleon sea is an- 
other nontrivial aspect of the nucleon structure. From 
the experimental side, some evidence has been found on 
a possible s — s asymmetry coming from global fits to 
data [5]. Notice that there is no fundamental symme- 
try preventing s{x) ^ s(x) in the nucleon, provided that 

/p [s{x) — s{x)] = 0. On the theoretical side, specula- 
tions about a possible \KH) component in the nucleon 
wave function, where K, H are virtual Kaon and Hy- 
peron states, leading to a s — s asymmetry, date since 
1987, with the pionering work by Signal and Thomas [7]. 
Since then on, several models have been proposed in the 
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literature [1111], with different predictions. It is worth to 
note also that a s — s asymmetry in the nucleon is gener- 
ated perturbatively starting at Next-to-Next to Leading 
Order (NNLO) TIT because at this order the splitting 
functions Pqq and Pqq are different. That asymmetry, al- 
though very small, should compete with the s — s asym- 
metry generated by the non-perturbative dynamics of the 
bound state. 

In this paper we shall investigate to what extent a non- 
perturbatively generated s — s asymmetry can describe 
the strange sea asymmetry found in recent global fits to 
experimental data. In order to do that, in Section |ll] 
we will revise a model for the generation of the intrin- 
sic s/s sea of the nucleon, following in Section 
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a comparison of the model with experimental data and 
recent parameterizations of the strange sea of the pro- 
ton. Section Hv] will be devoted to further discussion and 
conclusions. 
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II. THE NON-PERTURBATIVE STRANGE SEA 
OF THE PROTON 



FIG. 1: Sea q/q parton distribution in the proton as given by 
Eq.[2l 



We start by considering a simple picture of the nucleon 
in the infinite momentum frame as being formed by three 
dressed valence quarks - valons, v{x) - which carry all of 
its momentum [TTj . 



Vuix) 



/3(a„ + l,6„ + l) 
1 



x'^-'il-xp , (1) 



with X the fraction of momentum carried by the valon 
with respect to the proton. 

In the framework of the Meson Cloud Model (MCM), 
the nucleon can fluctuate to a meson-baryon bound state 
carrying zero net strangeness. As a first step in such 
a process, we may consider that each valon can emit 
a gluon which, before interacting, decays perturbatively 
into a ss pair. The probability of having such a perturba- 
tive qq pair can be computed in terms of Altarelli-Parisi 
splitting functions [2] 



4 1 + (l-z)2 
3 ~z ' 

I + (1 - ^r) 



(2) 



These functions have the physical interpretation as the 
probability of gluon emision and qq creation with momen- 
tum fraction z from a parent quark or gluon respectively. 
Hence, 



q{x,Ql)^N 
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[y^v(z) (3) 



is the joint probability density of obtaining a quark or 
anti-quark coming from subsequent decays w — )■ v + 5 
and g q + q a,t some fixed low Qg. In eq. ([3|, 
is a normalization constant which should scale with the 
masses of the flavors being created so that to a heavier 
flavor, corresponds a smaller N . Since the valon distribu- 
tion does not depend on Qq [IT , the scale dependence in 
eq. ([3| only exhibits through the strong coupling constant 
agf The range of values of Qq at which the process of 
virtual pair creation occurs in this approach is typically 
below 1 GeV^. A tipical sea quark distribution obtained 
from eq. ([3]), is shown in Fig. ([T]), where, in the spirit of 
keeping as simple as possible the description of the pro- 
cess , whe have assumed Vu = Vd and used — aa — Q-h, 
bu = bd = 2 in eqs 0. 

Once a ss pair is produced, it can rearrange itself 
with the remaining valons so as to form a most ener- 
getically favored Kaon-Hyperon bound state. In order 
to obtain the Kaon and Hyperon probability densities in 
the \KH) component of the proton wave function, the 
well known approach of the recombination model |12| 
has been used 8]. Notice however that the Kaon and 
Hyperon probability densities obtained in this way can 
be represented in terms of the simple forms [5] 



Pk{x) 



1 



-x"" {1 - x)'"' , 



l3{aK + l,bK + l) 

Ph{x) = ^ x--{l~xf", (4) 

l3{aH + l,bH + 1) 

which are both properly normalized to unity, as it should 
be for a bound state of one Kaon and one Hyperon and 
in order to cope with the zero net strangeness of the 
proton. These forms are also consistent with the valon 
model for a hadron made of two partons bound state. 
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The coefRcients ax, bx, clh, bn in eqs. Q are not 
independent. In fact, as the Kaon and Hyperon have to 
exhaust the momentum of the proton, then 



dx [xPh{x) + xPk{x)] = 1 



(5) 



giving the constraint 

r(aK + 



r(aK + l)r(aK + &K + 3) 
r(aH + + 2)r(aH + 2) 
r(a„ + l)r(a„ + 6h + 3) 



= 1 



(6) 



Finally, the non-perturbative strange and anti-strange 
sea distributions in the nucleon can be computed by 
means of the two-level convolution formulas 
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-NP 



(x) 



(x) 



N I ' ^Pniy) SH{x/y), 

y 

1 dy 



N 



y 



-PkIv) SK{x/y), 



(7) 

(8) 



where the sources sh{x) and sk{x) are the probabil- 
ity densities of the strange valence quark and anti-quark 
in the Hyperon and Kaon respectively, evaluated at the 
hadronic scale Qq |7j. In principle, to obtain the non- 
perturbative distributions given by eqs. ([S]), one should 
sum over all the strange Kaon-Hyperon fluctuations of 
the nucleon but, since such hadronic Fock states are nec- 
essarilly ofF-shell, the most likely configurations are those 
closest to the nucleon energy-shell, namely hPK'^ ^ S+if^ 
and TPK^ , for a proton state. For the sake of simplicity, 
we will only consider a generic Kaon and Hyperon inside 
the proton. 

For the sh{x) and sk{x) probability densities in eqs. 
([t]) and ([8| we also used the simple forms 



sh{x) 
sk{x) 



1 



/3(a,H + l,6,ff + 1) 
1 

P{asK + ^,bsK + 1) 



x"'" {I - xf'" , 
x''^" {1 ~ xf'" , (9) 



according to the valon model. The coefRcients 
o-sK, bsKi cisH, bsH in eqs. ([9| have to be determined 
by comparison to experimental data. 



III. COMPARISON WITH DATA 

In order to fix the coefficients of the model, we compare 
with the s and s parton distribution functions ^15J given 
in Ref . [6 , which have been recently obtained in a global 
fit to DIS data. The comparison has been done by means 
of a simultaneous least square fit of the model to xs{x) + 
xs{x) and xs{x) — xs{x) at = 20 GeV^. The results 
of the fit are shown in Figs. and ^ and in Table [U 

The procedure of the fit was as follows: i) a. set of 
initial values for the parameters was chosen and the non- 
perturbative s^^ and s^^ strange quark distribution 



TABLE I: Coefficients for s{x) and s{x) obtained from simul- 
taneous fits to the strange parton distributions funcions in 
the proton as given in Ref. [B]. uh = 2.889 is fixed by the 
requirement of momentum conservation given by eq. 



bn 


1.563 


±0.166 


clk 


2.403 


±0.127 


bK 


4.164 


±0.188 


a,sH 


0.669 


±0.052 


bsH 


8.539 X 10"^ 


±0.087 


dsK 


23.203 


±0.085 


bsK 


0.602 


±0.002 


N 


0.019 


0.001 




-0.004 



FIG. 2; xs{x) - xs{x) at = 20 GeV^. Dashed line: the 
asymmetry obtained in Ref. [6]. Full line: the result of our 
fit. The shadow area is the uncertainty of the fit. 



functions were calculated according to eqs. ([T]) and (|8|. 
Then ii) the s^^ and s^^ distributions were evolved 
from Qq up to — 20 GeV^ and the squared distance 



n 

E 

1=1 



{yd{x.)-yf{x,)y 

{y,{x,)-yf{x,)y 
alixi) 



was calculated. In eq. (10), yk{xi) 
-^^(Q2,x,) and yf 



xs 



NP 



(10) 



xs"" ((^"^Xi) and ?/^'° = xs{Q^,Xi) ± xs{Q'^,Xi), with 
index k — d, s refering to the difference and the sum 
respectively. An arbitrary error cr^/s(a:i) corresponding 
respectively to 10% of the value of the difference and the 
sum at Xi has been considered to perform the fit. Finally 
a new set of parameters was chosen and the procedure 
has been repeated until a minimum in 5^ was reached. 



4 




FIG. 3: xs{x) +xs{x) at = 20 GeV^ Dashed line: the 
distribution obtained in Ref. [6]. Full line: the result of our 
fit. In the insert is shown the result of our fit for 0.1 < x < 1. 
The uncertainty of the fit is represented by the shaded area. 
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FIG. 4: Full line: Hyperon probability density in the \KH} 
Fock state of the proton. Dashed line: Kaon probability den- 
sity in the \KH). 



The package MINUIT has been used to perform the fit. 
Notice also that the QCD evolution of the combination 
s + s depends on the whole set of parton distribution 
functions in the proton, for which we used the pdf's of 
Ref. 6^. 

The best fit has been obtained using Qq — 0.3 GeV 
as the starting point for QCD evolution. As can be seen 
in Figs. ([2| and ([s]), the model represents fairly well the 
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FIG. 5: Strange (left) and anti-strange (right) quark distribu- 
tions in the proton at = 20 GeV^. Our model compared to 
the V. Barone et al. MRST [S] and CTEQ5 [H] strange 
parton distribution functions. 



strange parton distributions found in Ref. [64 for x ^ 0.1, 
whereas for smaller values of x our results are below the 
results of the global fit of Ref. [B] . This can be due to a 
deficit in the content of gluons, as seems to be indicated 
by the fact that a good fit is only obtained for extremely 
low values of Qo- 

Concerning the model itself, the Kaon and Hyperon 
distributions functions in the \KH) Fock state of the pro- 
ton are displayed in Fig. (l4|. The momentum carried by 
the Hyperon in the \KH)Voc\ state is xPh{x) — 0.6 
while for the Kaon we obtained xPk{x) = 0.4, agree- 
ing with the common intuition that the Hyperon carries 
more momentum than the Kaon in the \KH) component 
of the proton wave-function. 

In Fig. ([5]) the strange and anti-strange quark distribu- 
tions at Q = 20 GeV^ are shown and compared to the 
xs and xs distributions found in Ref. [6j. We also com- 
pare to the MRST and CTEQ5 [II] strange quark 
pdfs, which have been determined imposing s = s. As 
shown in the figure, while our strange quark and anti- 
quark pdfs and those of Ref. [6] are consistent in the 
full range 0.01 < a; < 1, they deviate from the behavior 
shown by the MRST and CTEQ5 sets at x > 0.3. 



IV. CONCLUSIONS 

In this paper we have presented a model for the non- 
perturbative structure of the strange sea of the proton. In 
the model, the non-perturbative s^^ and s^^ intrinsic 
sea quark distributions of the proton are given in terms of 
a convolution of Hyperon and Kaon probability densities 
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and valence quark distributions inside a \KH) compo- 
nent of the proton wave function. This naturally gen- 
erates an asymmetry in the momentum distributions of 
the strange and anti-strange quarks in the proton, since 
the Hyperon, being heaviest than the Kaon, carries more 
momentum in the \KH) wave funtion component. The 
model depends on eight parameters which have to be 
fixed by fits to experimental data. 

Parameters of the model have been fixed by fits to the 
strange quark distributions found in a recent global fit to 
DIS data [6 . As shown in section III , the model describes 
qualitatively well the behavior of the s — s as well as s -I- s 
distributions, being this agreement better in the region 
X ^ 0.1. For lower a;, the s -\- s distribution is below the 
corresponding curve given by the distributions of Ref. 'B] . 
The probability of the \KH) fluctuation of the proton is 
about 0.02%, as given by the parameter N . 

The fact that the model does not describe well the s-l-s 
distribution at cc < 0.1 is expected by several reasons. 
First of all, the parton distributions found in Ref. [6 
were determined in a global fit to DIS data, instead our 
s and s distributions have been determined by fits to the 
s — s and s + s distributions found in [B] . This restrict 
the space of parameters allowed to the fit. Second, the 
gluon distribution found in Ref. [6] seems to be insuffi- 
cient to generate enough s/s quarks at low momentum, 
as evidenced by the fact that while the s — s distribution 
is well described in the whole range < a: < 1, the s-l-s 
is not. This could also explain why a good fit is obtained 
only for extremely low values of Qqj the scale at which 
perturbative QCD evolution starts. And third, a mean- 
ingful comparison of the model has to be done through a 



global fit to experimental data. 

Notice also that no NNLO effects have been included 
in the QCD evolution of the strange parton distribu- 
tions. The inclusion of those effects should produce a 
slightly bigger non-perturbative contribution to the pro- 
ton wave funtion to compensate the opposite asymmetry 
which arises at NNLO [TO] . 

Finally we would like to emphasize that the model pre- 
sented here can be taken as a phenomenologically moti- 
vated alternative for the description of the input strange 
and anti-strange quark distributions used in fits to ex- 
perimental data. Notice that the model allows for a full 
representation of the non-perturbative processes inside 
the proton in terms of well known mechanisms such as 
the splitting of quarks and gluons and recombination, 
which after all must account for the dynamics of the pro- 
ton as a bound state. Another interesting aspect of the 
model is that it can shed light on the parton structure 
of real Kaons and Hyperons, since the structure of the 
strange sea of the proton is related to the strange va- 
lence quark distributions at a low scale of the former 
strange mesons and baryons. 
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